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An experimental investigation was undertaken to measure the intrinsic elastic properties of several of 
the microstructural components of human vertebral trabecular bone and tibial cortical bone by the 
nanoindentation method. Specimens from two thoracic vertebrae (T-12) and two tibiae were obtained 
from frozen, unembalmed human male cadavers aged 57 and 61 years. After drying and mounting in 
epoxy resin, nanoindentation tests were conducted to measure Young's modulus and the hardness of 
individual trabeculae in the vertebrae and single osteons, and interstitial lamellae in the tibiae. 
Measurements on the vertebral trabeculae were made in the transverse direction, and the average 
Young's modulus was found to be 13.5 ± 2.0 GPa. The tibial specimens were tested in the longitudinal 
direction, yielding moduli of 22.5 ± 1.3 GPa for the osteons and 25.8 ± 0.7 GPa for the interstitial 
lamellae. Analysis of variance showed that the differences in the measured moduli are statistically 
significant. Hardness differences among the various microstructural components were also observed. 
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The elastic properties of bone have been measured at 
the microstructural level in many investigations 1 " 14 . 
Some of the microstructural features of interest 
include individual trabeculae 1 " 8 , single osteons 9 " 11 
and the thin cortical shell 12 " 14 . For many years, it was 
assumed that the elastic moduli of these microstruc- 
tui-al components are similar to those of dense cortical 
bone. However, recent investigations have suggested 
that this may not be the case; rather, the moduli of 
some of the components may be smaller 1-14 . In 
addition, it appears that significant variations in 
Young's modulus may exist within a given microstruc- 
tural component; for example, Young's moduli for 
individual trabeculae reported in the literature range 
from 1 to 14 GPa 18 . It should be noted, however, that 
even the largest of these values is still smaller than the 
frequently quoted modulus for dense cortical bone, 
E = 17.1 GPa 15 . 

Observations of specimen size influences on the 
measured modulus also suggest that the elastic moduli 
of the individual microstructural components of bone 
may not be the same as the macroscopically measured 
values. Choi et al. 2 found that Young's modulus for 
cortical bone measured in micro-bending tests 
decreases with decreasing specimen size. In separate 
studies, Choi et al. 2 and Lotz et ai." found that 
Young's modulus for cortical bone obtained from 
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micro-bending tests (5.4 and 12.5 GPa, respectively) is 
considerably smaller than that obtained by others in 
tensile tests of large specimens (17.1 GPa) 15 . One 
explanation for these observations is that the elastic 
properties of the microstructural components of bone 
are not the same as the macroscopic values. In this 
regard, a technique by which the elastic properties of 
the individual microstructural components can be 
measured would be of great value in understanding 
the mechanical behaviour of bone at the microstruc- 
tural level. 

Nanoindentation is one such technique. Developed 
over the last 10 years, nanoindentation is now used 
widely in the materials science community for probing 
the mechanical properties of thin films, small volumes 
and small microstructural features 16 . The properties 
most commonly measured are Young's modulus, E, 
and hardness, H 1617 . One of the great advantages of 
the technique is its ability to probe a surface and map 
its prnpertifis on a spatially-resolved basis often with a 
resolution of better than 1pm. Since many of the 
microstructural features of interest in bone are several 
micrometres or more in dimension, the nanoindenta- 
tion technique offers a means by which their intrinsic 
mechanical properties can be measured directly. In the 
present study, the nanoindentation method was used 
to determine the elastic modulus, and hardness of 
individual vertebral trabeculae and single osteons and 
interstitial lamellae in tibial cortical bone 



1325 



Biomaterials 1997, Vol. 18 No. 20 



1326 



Nanoindentation properties of bone: J.-Y. Rho et a/. 



MATERIALS AND METHODS 

Two thoracic vertebrae (T-12) and two tibiae were 
obtained from two frozen, unembalmed human male 
cadavers 57 and 61 years of age (one vertebra and one 
tibia from each cadaver]. For the vertebrae, the 
spinous processes were first removed, and the 
vertebral body was cut to a thickness of 2 mm along 
the sagittal plane with a low-speed diamond saw 
(Buehler Isomet 1000, Buehler, Lake Bluff, IL, USA) 
under constant deionized water irrigation. Each 
section was obtained from the anterior aspect. 
Following removal of the bone marrow using a water 
jet (Teledyne Water Pic, Fort Collins, CO, USA), the 
sections were dehydrated in a series of alcohols to 
avoid problems that can be caused by surface liquid 
films in identifying the point of first contact during 
indenter approach to the specimen surface. With 
machine and software modifications currently under 
development, such problems can potentially be 
overcome in future work, and it may be possible to 
test wet as well as dry specimens. Influences of 
specimen dehydration on nanoindentation property 
measurement will be discussed in the next section. 

After dehydration, the specimens were embedded 
in epoxy resin (Epo-Thin M low viscosity epoxy; 
Buehler) to provide support for the porous network 
and then metallographically polished to produce the 
smooth surfaces needed for nanoindentation testing, 
first with silicon carbide abrasive papers of 
decreasing grit size (600, 800 and 1200 grit), and 
finally with diamond suspensions (0.3 and 0.05 /mi 
particle size) embedded in soft polishing cloths. The 
specimens were washed in deionized water between 
each polishing step to remove debris, taking care that 
the surface of the bone was not demineralized. The 
polished surfaces of the vertebral sections were 
oriented such that nanoindentation testing measured 
their properties in the transverse direction, 
Longitudinal properties were not measured due tu 
difficulties in finding and unambiguously identifying 
the individual lamellae in the polished sections. 
Sections of the tibiae were prepared in a similar 
fashion, but with the specimen oriented for 
measurements in the longitudinal direction, for 
which individual osteons and trabeculae were easy 
to identify. 

All experiments were performed using the Nano 
Indenter II (Nano Instruments, Inc., Knoxville, TN, 
USA) at the Oak Ridge National Laboratory. This fully 
automated hardness testing system makes small 
indentations at precise positions on a surface while 
< f i i < fik 1 monitorin *>> In id and displacements 
of the indenter with resolutions of 0.3 fiN and 0.16 nm, 
respectively. A sharp Berkovich diamond indenter, a 
three-sided pyramid with the same area-to-depth ratio 
as the Vickers indenter used commonly in conventional 
microhardness testing, was used for all measurements. 
In a typical experiment, the microstructural feature to 
be examined was located in the microscope and 
positioned beneath the indenter using the x-y table. 
The indenter was then slowly driven toward the 
surface at a constant displacement rate of lOnms" 1 
until surface contact was detected by changes in the 
load and displacement signals. After contact, a 
permanent hardness impression was made by driving 



the indenter into the specimen to a depth of 1000 nm 
at a constant loading rate of 750 /iN s 1 , holding at this 
load for a period of 10 s and then unloading to 15% of 
the peak load at a rate equal to half that used during 
loading. At the end of the unloading cycle, the indenter 
was held on the surface for a period of 100 s to establish 
the rate of thermal drift in the machine and specimen 
for correction of the data, and then completely 
wit! i n \ i oximatel) £h indentations were 
made in each feature. Any indentations close to the 
mounting resin were removed from the data set to 
minimize effects of embedding on the measurements. 

The indentation load-displacement data obtained in 
these tests were analysed to determine the hardness, 
H, and Young's modulus, E, using the method of 
Oliver and Pharr 1 ". The first step in the analysis 
procedure is to determine the contact stiffness, S, a 
measure of the resistance of the bone to elastic 
deformation. As shown in the schematic depiction of a 
typical indentation load-displacement curve during 
one cycle of loading and unloading in Figure 1, the 
contact stiffness is experimentally measured from the 
load-displacement data as the slope of the upper 
portion of the unloading curve. Young's modulus is 
related to the contact stiffness by 

where Eb and v b are Young's modulus and Poisson's 
ratio, respectively, for the bone, Ei and Vj are the same 
quantities for the indenter, and the factor ft is a constant 
which depends on the indenter geometry. For the 
Berkovich indenter, 0= 1.034. This equation is 
derived under the assumption that the material is 
homogeneous and isotropic in its elastic properties. 
When used to measure the elastic modulus of an 
anisotropic material, the modulus derived from 
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Figure 1 A schematic illustration of typical indentation 
load-displacement behaviour during one cycle of loading 
and unloading, showing quantities used to determine 
hardness and Young's modulus from the data. The quanti- 
ties are the peak indentation load, P max , the maximum 
indenter displacement, h max , and the initial unloading 
stiffness, S. 
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Equation 1 is an average of the anisotropic elastic 
constants which is biased toward the modulus in the 
direction of testing 19,20 . If the elastic symmetry of the 
material is known, it is possible to perform experiments 
to determine the individual anisotropic elastic 
constants. However, in this paper, only the simple 
average value is reported. For a diamond indenter, 
Ei = 1141 GPa and Vj = 0.07, and we assume that 
Poissun's ratio for bone is v h = 0.3 1 '\ A sensitivity 
study showed that varying v b in the range 0.2-0.4 
changed the measured values of by no more than 
8%. 

The use of Equation 1 in the measurement of elastic 
modulus requires an independent measure of the 
contact area, A. As detailed elsewhere 18 , the contact 
area can be estimated directly from the indentation 
load-displacement data through a process which 
requires a precise knowledge of the shape of the 
indenter determined in special calibration experiments 
and measurement of the 'plastic depth' of the contact, 
also referred to as the 'contact depth'. The contact 
depth can be computed from the contact stiffness, S, 
the peak load, P mal , and the maximum depth of 
penetration, £ max , all of which are experimentally 
measurable from the load-displacement data. The 
hardness, H, which physically represents the mean 
pressure the material can support, can also be 
determined from this procedure. The hardness is 
defined as the maximum load, P maK , divided by the 
projected area of the contact impression, A, i.e. 

H = PU (2) 

Determining the contact area in this way has the 
advantage that both E and H can be measured in one 
simple test. In addition, mechanical property 
measurements can be made without having to imagine 
the indentation. In metals and ceramics, it is not 
uncommon for E and H to be measured using 
indentations only tens of nanometres deep 18 . 

A total of 199 indentations were produced in the 
study. In the vertebrae, three to eight indentations 
were made in 14 separate trabeculae. Fur (he tibiae, 
12 osteons and six interstitial lamellae were tested 
using three to five indentations in each osteon and 
five to 10 indentations in each interstitial lamella. 
Significant differences in the elastic properties of 
trabeculae, osteons and interstitial lamellae were 
analysed using one-way analysis of variance 
(ANOVA). Scheffe's test was then employed to find 
differences between the elastic properties of the 
vertebral trabeculae and dense tibial cortical bone 
(osteons and interstitial lamellae). 

A fused silica calibration standard was tested to 
check on the quality of the measurements. Fused silica 
is often used for calibration in nanoindentation testing 
because its relatively low modulus-to-hardness ratio, 
E/H, leads to a large amount of elastic recovery during 
unloading, and this improves measurement accuracy 18 . 
The elastic modulus of the fused silica sample was 
measured to be 71.8 ±0.4 GPa, which is within 0.3% 
of the known value, 72.0 GPa. This provides some 
assurance that the machine was properly calibrated 
and the analysis procedures were properly 
implemented. 



RESULTS 

Figure 2 is an optical micrograph of a linear array of 
nanoindentations made in one of the tibiae, illustrating 
how both osteons and interstitial lamellae can be 
individually sampled. The indentations shown in the 
figure are approximately 1 ^m deep and 5 ^m on an 
edge. A typical set of load-displacement data for an 
osteon is shown in Figure 3. The general shapes of the 
loading and unloading curves are similar to those 
observed in metals and ceramics 18 . 

A summary of the elastic moduli and hardness for the 
various microstruntiiral components of the vertebrae 
and tibiae is presented in Table 1. All data in the table 
were obtained at an indentation depth of about 
1000 nm (l^tm), It is important to remember that 
measurements on the trabecula were made in the 
transverse direction, while those on the cortical bone 
were in the longitudinal direction. Also included in 
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Figure 2 An optical micrograph of a linear array of nanoin- 
dentations in tibia! cortical bone. Some indentations are 
located in an osteon while others are in an adjacent intersti- 
tial lamella. 
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Figure 3 Load-displacement data from a nanoindentation 
test in a tibial osteon. 
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Table 1 Elastic moduli and hardnesses for the various bone microstructural components 
Bone type; feature Specimen 
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the table are the values of E and H for each microstruc- 
tural component averaged over all indentations in both 
specimens. The average results show that Young's 
modulus is greatest for the tibial interstitial lamellae 
(is = 25.8 GPa), followed in decreasing order by the 
tibial osteons [E = 22.5 GPa) and the vertebral 
trabeculae [E = 13.4 GPa]. The hardnesses follow in 
the same order; the hardest component is the tibial 
interstitial lamellae (H == 736 MPa), followed by the 
tibial osteons (H = 614 MPa) and the vertebral 
trabeculae (H = 468 MPa I. ANOVAs showed that the 
mean values of E and H are statistically different 
(P < 0.0001) for all the bone components. 



DISCUSSION 

The nanoindentation technique has several advantages 
over conventional microhardness measurement 
methods, such as the Vickers and Knoop techniques. 
The major benefits accrue from the accurate positioning 
capability (better than l,um) and the high resolution 
load and depth-sensing capabilities which enable 
relatively small areas of material to be tested and 
measured. Also, unlike conventional microhardness 
testing, it is not necessary for the area of nanoindenta- 
tions to be optically measured in order for the hardness 



or elastic modulus to be determined. In this way, 
indentations smaller than the limit of optical resolution 
can be used in the testing 21 , and very small features can 
be explored. In this study, osteons, interstitial lamellae 
and individual trabeculae with characteristic 
dimensions of tens to hundreds of micrometres were 
examined 22,23 . 

The average Young's modulus for vertebral 
trabeculae measured in this study, £" = 13.4 GPa, is 
considerably smaller than that for the components of 
dense cortical bone, i.e. F = 22.4 GPa for osteons and 
E = 25.8 GPa for interstitial lamellae. However, some 
caution should be exercised in interpreting this result, 
since Young's moduli for the trabeculae were measured 
in the transverse direction, while those for the cortical 
bone (tibial osteons and interstitial lamellae) were 
measured in the longitudinal direction. In this context, 
it is notable that Roy et al. 2 *, using nanoindentation 
methods similar to those employed here, recently 
characterized both the transverse and the longitudinal 
moduli of single trabeculae and found a significant 
elastic anisotropy. The transverse modulus measured 
for horizontally-oriented trabeculae, 16.0 GPa, is 
similar to that measured here. However, the modulus 
they found in the longitudinal direction 22.7 GPa, is 
significantly higher, and in fact is very similar to the 
current measurements for tibial osteons. Thus, when 
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elastic anisotropy is taken into account, the differences 
in moduli between trabecular and cortical specimens 
measured in this study may not be as significant as 
they first appear. 

In combination with the results of Roy et alS*, the 
nanoindentation results obtained here suggest that the 
moduli of individual trabeculae are significantly 
greater than the widely used value, E = 5.4 GPa, 
obtained by microbending testing 2 . Some of this 
discrepancy could be due to the drying of the 
specimens prior to the nanoindentation testing. To 
estimate how large this effect might be, it is useful to 
refer to the work of Townsend and Rose 8 , who found 
that drying increases Young's moduli of individual 
trabeculae by about 24%. Assuming this same factor 
applies to the nanoindentation measurements, the 
transverse modulus for individual trabeculae measured 
in this study would be about 10.2 GPa after 
approximate correction for the effects of drying, and 
the longitudinal modulus reported by Roy et al. 2 * 
would correct to about 18.3 GPa. Thus, even when 
drying is taken into account, the nanoindentation 
results still suggest a higher modulus for trabeculae 
than what is usually assumed based on microbend 
testing. The reason for this discrepancy is not clear, 
but it could arise from difficulties encountered in 
making accurate mechanical property measurements 
with small bend specimens. These problems include: 
(1) influences of microstructural defects such as 
cement lines and voids (Haversian and Volkmann 
canals, lacuna, osteocytes. canaliculi and so on] on the 
measured displacements: (2) uncertainties in specimen 
geometry, which are often exacerbated at small scales; 
and (3) problems in properly seating and aligning 
small bend specimens in small test fixtures. These 
problems are largely avoided by the nanoindentation 
method, since testing can be done in areas specifically 
chosen to avoid defects, the test geometry is determined 
by the shape of an accurately calibrated diamond 
indenter. and seating problems are minimized through 
well-developed surface contact detection schemes 
built into the software that control the nanoindentation 
testing system. It is thus suggested that the inherent 
modulus of individual trabeculae may be considerably 
higher than is commonly believed. 

Another question which arises is the extent to which 
embedding in epoxy affects the nanoindentation 
measurements. To minimize these effects, indentations 
close to the bone-resin interface were not included in 
the data set. Some insight into this problem can be 
gained from the work of Evans et a/. 25 , in which a 4% 
difference in hardness between embedded and 
unembedded bone was observed. However, as noted 
by the authors, it was not possible to separate the true 
effects of embedding from those caused by dehydration 
during the embedding process. Further studies of the 
effects of dehydration and specimen embedding on 
nanoindentation proj rt m isurement are warranted. 
For cortical bone, it should be feasible to make 
nanoindentation measurements on both wet and dry 
bone in both embedded and unembedded specimens 
and thus directly address both influences at once. On 
the other hand, it seems unlikely that trabecular bone 
can be studied without embedding, so answers to these 
questions may have to come from tests on cortical bone. 

To the best of our knowledge, the elastic moduli of 



interstitial lamellae have not been measured 
previously. The results obtained here show that there 
is a statistically significant difference in Young's 
moduli of interstitial lamellae and osteons in the same 
bone, with the modulus of the interstitial lamellae 
being about 15% greater. Based on these observations, 
the modulus of a macroscopic sample of cortical bone 
should fall somewhere between that of the osteons and 
interstitial lamellae, with tho exact value depending on 
the microstructural arrangement of the two components 
and their relative volume fractions. If the volume 
fractions were known, isostress or isostrain averages 
could be used to provide bounds on the macroscopic 
modulus. The isostrain average would probably give 
the better estimate, since the osteons and interstitial 
lamellae are arranged in the longitudinal direction in a 
manner which subjects them more to isostrain than 
is tfress conditions As to -why Young's moduli of 
osteons and interstitial lamellae are different, it has 
been observed using backscattering scanning electron 
microscopy that interstitial bone is more mineralized 
than osteons 26 . The lower elastic modulus for osteons 
may thus result from the osteons being newer bone 
material than the interstitial lamellae, since newer 
bone is known to have lower mineral content 27 . As 
further support for this notion, Evans and Vincentelli 28 
have reported that osteons tend to reduce both the 
strength and the elastic modulus of cortical bone. 

This initial study demonstrates that the nanoindenta- 
tion method can offer valuable insight into tho elastic 
properties of the microstructural components of bone. 
Properties measured by nanoindentation could prove 
useful in the development of theoretical micromecha- 
nical models 29 ^ 31 and in finite element modelling. The 
nanoindentation characterization of the mechanical 
properties of bone at the microstructural level could 
also provide information useful in understanding the 
complex interactions of biological and cellular 
mechanisms for osseous disorders. 
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